jth p erm i Symposium : Monterey, CA : 28 Oct- 2 Nov 2012 



o 

(N 



in 

(N 

X 

Ok 

6 

CO 

(N 
> 

OS 

o 
in 

cn 
o 
m 



Constraints on the Intergalactic Magnetic Field from Gamma-Ray 
Observations of Blazars 

J. D. Finke 

Space Science Division Code 7653, US Naval Research Laboratory, 4555 Overlook Ave. SW, 20375, 
USA 

L. Reyes 

Physics Department, California Polytechnic State University, San Luis Obispo, CA 94307, USA 
M. Georganopoulos 

Department of Physics, University of Maryland Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250, 
USA and 

NASA Goddard Space Flight Center, Code 660, Greenbelt, MD 20771, USA 
on behalf of the Fermi-LAT Collaboration 



Gamma rays from distant blazars interact with the extragalactic background light, creating electron-positron 
pairs, and reducing the gamma-ray flux measured by ground-based atmospheric Cherenkov gamma-ray tele- 
scopes. These pairs can Compton-scatter the cosmic microwave background, creating a gamma-ray signature 
observable by the Fermi Large Area Telesope (LAT) . The signature is also dependent on the intergalactic mag- 
netic field (IGMF), since it can deflect the pairs from our line of sight, reducing the gamma-ray emission. We 
present preliminary constraints on the IGMF using Fermi-LAT and Cherenkov telescope observations, ruling 
out both very large and very small values of the IGMF strength. 



1. Introduction 

The extragalactic backgroung light (EBL) from the 
infrared (IR), through the optical and into the ul- 
traviolet (UV) is dominated by emission, either di- 
rectly or through dust absorption and reradiation, of 
all the stars which have ever existed in the observ- 
able universe. It contains information about the cos- 
mological expansion, star formation history, dust ex- 
tinction and radiation in the universe, and so can 
potentially provide constraints on a number of cos- 
mologically interesting parameters. However, its di- 
rect detection is hampered by the bright foreground 
emission from the Earth's atmosphere, the solar sys- 
tem, and the Galaxy. This can be avoided by mea- 
surements fro m spacecraft outside the Earth's at- 
mosphere [e.g., Bernstein 2007 . Bernstein et all 2002 . 
iDwek et al l Il998l lHauser et all Il998l iMattilaf booal 
or solar system Edclstei n et all l2000l iLeinert et al 
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ci)(M-.ries (MeV-GeV) by EGRET [Hartman et all 
and another at very-high energies (VHE; 



> 0.1 TeV) by ground-b ased atmospheric C herenkov 
telescopes [ACTs; e.g., iPunch et aflll992j . It was 
almost immediately realized that 7-ray observations 
of extragal actic blazar s could be used to constrain 
the EBL iBiller et all Il995l IDwek fc Slavinl 11994 . 

cerl 



Stecker et all Il992j . 



Madau fc Phinnevl Il996l IStecker fc de Jage 



1993, 



The detection of the hard 
the BL Lac 1ES 1101-232 with 
put strong constraints on the EBL, 
that predict high opacity if one 
VHE spectrum cannot be 
differential pho- 



1998, lMatsuoka et al.ll201lllMurthv et alfeOOlllToUer 



spectrum from 
HESS seems to 
ruling out models 
assumes the intrinsic 
harder than T = 1.5, where the 
ton f lux $(£■) = dN/dE oc E~ T [Aharonian et all 
|2006|. However, the T = 1.5 constraint has been 
questioned, and a number of theoretical possibilities 
have been ra ised that could account for harder VHE 
spect ra lAharonian. Khangulvan. fc Costamantel 



1983], or through galaxy cou nts, which in genera l ; 
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Bottcher. Dermer, fc FinkeT 
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Stecker. Baring, fc Summerlinl 2007| . Nonethe- 



lower limits to the E BL [e .g-.lMad au & Pozzetti pOO 1 
iMarsden et aTll2009j . See lHauser fc Dwekl (200?| for a 
review of EBL measurements, constraints, and mod- 
els. 

Shortly after the discovery of the cosmic microwave 
backg round (CMB) radiation Penzias fc Wilsonl 
Il965| . it was realized that this and other radiation 
fields would interact with extragalactic 7-rays, pro- 
ducing e + e~ pairs, ef f ectively absorbing the 7 rays 
Fazio fc Stecker! I1970L iGould fc Schreden Il967| . In 



less, the basic idea, that the intrinsic VHE photon 
index cannot be harder than a certain number, 
has been used by many authors to constrain the 



EBL [e.g.. lAharonian et all l2007cl lAlbert et all 



the 1990s, extragalactic high-energy 7-ray astronomy 
took two major leaps forward, one with the launch 
of the Compton Gamma-Ray Observatory and the 
first detections of extragalactic sources at high 



20081 IFinke fc Razzaaud 120091 (iVIazin fc Rauej 120071 
Schroedterll2005| . 

The new era of 7-ray astronomy, began with the 
launch of the Fermi Gamma-Ray Space Telescope has 
brought additional constraints on the EBL. At the 
energies observed by the LAT, the extragalactic 7- 
ray sky below ~ 10 GeV is expected to be entirely 
transparent to 7-rays, at least back to the era o f re- 
comb ination [z ~ 1000; e.g.. IFinke et all l2010l lOhl 
|2001| . while photons observed in the 10 GeV - 300 
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GeV range should be attenuated by UV/optical pho- 
tons if they originate from sources at z > 0.5. This 
suggests a possible way to constrain the EBL. The 
LAT spectrum below 10 GeV can be extrapolated to 
higher energies, and should be an upper limit on the 
intrinsic spectrum in the range where the EBL at- 
tenuates the 7-rays. That is, everywhere for the 7- 
ray spectrum d 2 log($)/(dlog(i?)) 2 < (the spectrum 
is concave). This seems to be a reasonable assump- 
tion, since no convex 7-ray spectrum has been ob- 
served, and it is difficult to imagine theoretical ways 
to pr oduce one, a l thoug h see Section [5] for a discus- 
sion. lAbdo et al.l 2010j | have used this technique to 
put upper limits on the optical/UV EBL absorption 
optical depth (t 77 ) for sources z > 0.5, and rule out 
with high significance some models which predict high 
t 77 with data from the first 11 months of LAT opera- 
tion from 5 bla zars and 2 gamma-ray bu rsts (GRBs). 
More recently, lAckermann et al" I |2012aj used a sim- 
ilar technique in a composite fit to 150 BL Lacs 
with 46 months of LAT data to constrain the high-z 
EBL even fu rther, and found agreement with most re- 
cent models iDommguez et al.l2011 , Finke et alj2010 , 



teractions Compton-scatter the CMB photons, pro- 
ducing GeV 7-ray emission which could themselves 
be absorbed by i nteractions with the EBL, pro- 



Gilmore et all 120091 



2012 



Franccschini ct al. 2008 
Kneiske k Doldl2010j . 

Below 100 GeV, the universe is expected to 
be transparent out to z ~ 0.1, although VHE 
photons from this rcdshift range should be at- 
tenuated by interactions with IR E BL photons. 
iGeorganopoulos. Finke. k Revesl |2010| sug gested a 



very similar technique to Abdo et al. (2010j |. applied 
to the VHE range. For those sources detected by both 
LAT and an ACT, one can extrapolate the LAT spec- 
trum (which should be unattenuated for these sources 
at low z) into the VHE range, and use it as an up- 
per limit on the intrinsic VHE flux. As with the 
LAT-only case, a comparison of this upper limit with 
the observed VHE flux allows one to compute a n up- 
per limit on r 77 . IGeorganopoulos et al.l [2010j used 
this technique to show that models which predicted 
high t 77 in the VHE range were strongly disfavored. 
Another possible way to estimate the intrinsic spec- 
trum of a source, and thus estimate r 77 , comes from 
modeling the full radio to GeV 7-ray spectral energy 
distribution (SED) of a 7-ray blazar with a standard 
synchrotron/synchrotron sclf-Compton (SSC) model. 
The SSC spectrum can be extrapolated to the VHE 
regi me, and compared with observations to estimate 



'77 



Mankuzhivil et alj201l| . [Pominguez et all |2013 [ 
has applied this technique to a sample of ~ 15 LAT 
and ACT-detected blazars to constrain the cosmic 
7-ray horizon, i. e., the energy where tv, = 1 for a 
certain redshift. iDwek k Krennrichl [2012] present a 
comprehensive review of recent attempts to constrain 
the EBL with 7-ray observations. 

These constraints come with caveats, however. 
One possibility is that the electron positron pairs 
that are produced by the 7-ray-EBL photon in- 



ducing a ca s cade [Aharonian et al. 1 11994 iDai et al.l 
l2002llPlagal 1 19951 ?]. If the intergalactic mag- 
netic field (IGMF) strength is low, the pairs will 
not be significantly deflected from our line of sight, 
and this could prod uce an observable feature i n the 
LAT bandpass [e.g., iNeronov k Semikod 12009] . In- 
deed, recently several authors have even used the 
non-detection of these cascades to put lower lim- 
its on the IGMF str ength [e.g.. iDolag et al.l l201l[ 
INeronov k Vovldl20ld . lTavecchio et alj|2010j . In gen- 
eral these efforts depend on the fact that emission 
from TeV blazars is relatively constant over long pe- 
riods of time. TeV va riability has not been observed 
in some blazars, [e.g.. lAharonian et al. I I2006L l2007cl ] ■ 
although ma ny blazars are highly v ariable at these 
energies [e.g.. lAharonian et al. 2007 a] . Observing and 
studying the time-dependent EBL-induced pair cas- 
cades from GRBs and blazar flares has be en suggested 
as a way to probe the IGMF parameters llchiki et al.l 
l2008l iMurase et al.ll2008l iRazzaaue et alj|2004| . The 
possibility of variable TeV emission has led to caveats 
in interpreting the IGMF cons traints from appar- 
ently non-variable TeV blazars [Dermer et al.l l201ll 
iTavlor et al.l 12011 1 . Other uncertainties such as the 
EBL intensity and TeV spectra errors, blazar jet ge- 
ometry and Doppler facto r, can further dec rease the 
lower limits on the IGMF [Arlen et al.1l2012 |. 

In this paper we report on our efforts to constrain 
both the EBL absorption optical depth, t 77 and the 
IGMF strength (Bjq) and coherence length (Xb) us- 
ing 7-ray observations of the blazars 1ES 1101—232 
and 1ES 0229+200 from both LAT and ACTs. This 
tec hnique can be seen as an extension of pre vious work 
by IGeorganopoulos. Finke. k Revesl |2010l |. We make 
use of data from the first 3.5 years of LAT operation, 
the analysis of which is described in Section [21 and 
ACT spectra from the literature. Our technique for 
constraining the IGMF is described in Section [3] We 
report our preliminary constrains on the IGMF pa- 
rameters in Section[3J Finally we discuss these results, 
their caveats and implications in Section [5] 



2. LAT Analysis 

To determine the LAT spectra of 1ES 1101-232 and 
1ES 0229+200 we considered all LAT events collected 
since the start of the mission for 3.5 years of operation 
(i.e. until 2012 February 16), a significant increase in 
statistics with respect to previous efforts. The data 
were analyzed using an official release of the Fermi 
ScienceTools (v9r27) and Pass 7 instrument response 
functions. Only photons satisfying the CLEAN event 
selection are included in the analysis in order to re- 
duce the likelihood that charged-particle background 
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particles are mistaken as high-energy photons. As 
usual, the spectral analysis of each source is based on 
the maximum likelihood technique using the standard 
likelihood analysis software. In both cases we con- 
sidered all the point-like sources within 15° from the 
source position (as determined from the 2FGL catalog 
and improved diffuse background components (Galac- 
tic and extragalactic) . A likelihood ratio test was used 
to find the best spectral model (power-law, broken 
power law, and log parabola) that fits the data. In 
both cases, there was insufficient evidence for devia- 
tion from a power-law spectrum. 

In addition to the standard maximum likelihood 
analysis performed to find the best fit to the data, 
the Log-of-the-likelihood (LL) profile as a function of 
the source's flux normalization (Fq) was calculated in 
order to fully characterize the uncertainty on this pa- 
rameter. Thus, instead of assuming a perfect Gau- 
sian distribution for the error of Fq and using 1, 2, 
and 3 standard deviations to calculate the 68%, 95% 
and 99% confidence intervals, we used the LL profile 
to calculate the actual confidence intervals assuming 
that — 2A(LL) is distributed as the chi-square prob- 
ability distribution with one degree of freedom. We 
found that this approach is necessary in order to cor- 
rectly determine the flux probability distribution of 
weak LAT sources such as IES 0229+200 and 1ES 
1101-232. 



3. Method for Constraining Models 

Our tech nique for constraining the E BL and IGMF 
is based on iGeorganopoulos et al.l [2010j j . with exten- 
tions including a sophisticated Monte Carlo (MC) 
technique to accurately determine the significance of 
the constraints, and the addition of a pair cascade 
component. 



3.1. Technique Assumptions 

We make the following assumptions, the first 
three of which are idential to the ones of 
IGeorganopoulos et al.l 2010j |: 

#1 We assume the MeV-TeV flux from BL Lacs in 
our sample are produced cospatially from the sources 
themselves (and not UHECR interactions in inter- 
galactic space), and that they are never convex (i.e., 
never is d 2 \og(<P) / (dlog(E)) 2 > 0). 

#2 We assume the objects are not variable at 7- 
rays within the statistical uncertainties of the mea- 
surements. Indeed, we have selected sources for our 
sample for which 7-ray variability has not been re- 
ported, in either the LAT or ACTs. 

#3 We assume the 7 rays will not avoid being ab- 
sorbed on their path to Earth, by converting to axions 



or some other exotic mecha nism [ e.g. Jde Angelis et al.l 
120071 ISanchez-Conde et aDl2009| 

#4 We assume pairs created by 7-ray-EBL pho- 
ton interactions will Compton-scatter the CMB, and 
will lose energy primarily through scattering and 
not through intergalactic plasma beam instab ili- 
ties [Broderick et~al]|2012llSchlickeiser et al.ll2012afbl |. 
Our technique for includin g the cascade comp onent is 
similar to the one used by I Meyer et al.l |2012j . 

We discuss the viability of these assumptions in Sec- 
tion [5] 



3.2. Ruling out Models 

Our technique for ruling out a model is illustrated 
in Figure [T] and has the following steps: 

Step 1. Select the model we wish to rule out. This 
includes selecting an EBL model (taken from the liter- 
ature), selecting an IGMF strength (Bio) an d coher- 
ence length (As); and selecting a blazar jet radiation 
opening angle and a "blazar lifetime" (tuazar), that is, 
the length of time the blaza r has been emitting 7-rays 
with its current luminosity [Dermer et al.ll201l| . 

Step 2. Given the integrated LAT flux (Flat) and 
photon index (r) and their errors for a particular 
blazar (see Section [5]), we draw a random Flat and T 
from a probability distribution function (PDF) which 
represents their errors. 

Step 3. For each energy bin of the VHE spec- 
trum, with a measured flux and error, we draw a ran- 
dom flux, Fvhe, assuming the flux errors are dis- 
tributed as a normal distribution. Each randomly 
drawn Fvhe(E) is deabsorbed with the EBL model 
we are testing to give an intrinsic flux, FvHE,int(E) = 
exp(r 77 (£')) x F VHE (E). 

Step 4- From FyHE,int(E), the contribution of the 
e + e~ pairs Compton-scattering the CMB, F, 



cascade 



is calculated. The deabsorbed emission is assumed 
not to extend beyond the highest energy bin, the 
most conservative assumption we can make. The cas- 
cade flux F ca S c a de is calc ulated u s ing th e formula of 
iDermer et ail [201 lj and iDermerl |2012>. The minor 
corrections to this formula from iMever et al.l [2012} 
should have no real effect on the results, since they 
affect the lowest flux portion of the cascade spectrum. 

Step 5. The randomly drawn LAT power-law spec- 
trum from step 2 is extrapolated to the VHE regime. 
For this MC iteration, the model is considered rejected 
if one of two criteria are met: (i) the cascade flux 
Fcascade exceeds the randomly drawn LAT integrated 
flux Flat', or (ii) any one of the deabsorbed flux bins 
from FvHE.int(E) exceed the extrapolated LAT flux, 
FLAT.ext, unless O.OIFlat < F cascade < Flat, in 
which case the model is never rejected for this iter- 
ation. If the cascade flux makes up a significant frac- 
tion of the observed LAT flux, we do not believe it 
can be extrapolated to the VHE regime and used to 
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PRELIMINARY 




E [TeV] 



Figure 1: This figure illustrates many of the steps in our 
method for ruling out models from Section |3. 21 using the 
7-ray spectrum for 1ES 1101-232. The LAT spectrum is 
shown as the bowtie, along with this spectrum 
extrapolated to the VHE regime as the dashed curves. 
The observed HESS spectrum is shown as the filled 
squares. The randomly drawn HESS points shown as 
empty diamonds (Fvhe) and the randomly drawn LAT 
spectrum is shown as a line, both of which are labeled 
"Step 2". The deabsorbed points are shown as the circles 
(FvHE,int) and labeled "Step 3". The cascaded 
component and the interpolated VHE spectrum used to 
calculate it are shown as dashed curves labeled "Step 4". 
The LAT spectrum extrapolated into the VHE regime is 
shown as the dashed line labeled "Step 5". For the MC 
iteration shown here, the model is ruled out by both 
criteria in Step 5, since Flat < F casca d e and for several 

points F LA T,ext < F V HE,int- 



constrain models. The rejection criteria in th i s step 
are based on ones from iGeorganopoulos et al.l [2010] 
and lMever et all |2012j . 

Step 6. Steps 2-5 above are repeated N tr iais times 
(we use Ntriais = 10 6 ) and the number of times the 
model is rejected N re j ect is counted. The probability 
the model is ruled out is P re ject = N reject /N tna i s . 



o. 
-6[ 
-8 

1 " 1 ° 
I "12 

a -14 
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Figure 2: The values of parameter space of B and Lb 
which are ruled out by our method, using LAT and 
HESS observations of 1ES 0229+200. The significance 
that the values are ruled out are given by the colors, in 
number of sigma, as indicated by the color bar. 



and 1ES 1101-232 can be found in Figs. [2] and El 
respectively. The probabilities that values of parame- 
ter space are ruled are converted to number of sigma, 
assuming the errors are distributed normally. High 
values of Bjq are ruled out by the constraint that de- 
absorbed TeV flux cannot exceed the LAT spectrum 
spectrum extrapolated into this regime; and low val- 
ues of B are ruled out by the constraint that the cas- 
cade cannot exceed LAT flux (see Section EOl step 5). 
A black strip in these figures indicates allowed values 
of Big an d Lb- Here the cascade flux was less than 
the observed LAT flux, but > 0.01 of the LAT flux, so 
that the LAT flux could not be simply extrapolated 
into the TeV regime. 

The combined results for both 1ES 0229+200 and 
1ES 1101-232 are shown in Fig. gj High values of 
Bjg are ruled out at greater than 5c, while lower 
values are ruled out at greater than 3a. 



4. Constraints on IGMF 

We apply the method outlined in Section l3~2l to 7- 
ray observations of 1ES 0229+200 and 1ES 1101-232. 
We use the HESS spe ctra for these sources from 
lAharonian et al.l [2007d | and lAharonian et al.l [2007bj |. 
respectively, and the LAT spectra calculated as out- 
lined in Section p We used the EBL model from 
iFinke et al.l [201o| 7a jet radiation opening angle of 0.1 
rad, and a blazar lifetime of tuazar = Hq , where Hq 
is the Hubble constant; that is, we assume the blazar 
has been emitting for essentially the entire age of the 
universe. Neither of these objects have been shown to 
be variable at 7-ray energies, with either the LAT or 
TeV instruments. 

The results from the application to 1ES 0229+200 



5. Discussion 

We have used LAT and TeV data to constrain 
IGMF parameters. Our preliminary results indicate 
high values of Big are ruled out with > 5a signifi- 
cance. This means that either one of our model as- 
sumptions is not correct f Section 13. 1[) . one of our pa- 
rameters is not correct, or the cascade flux must con- 
tribute at least 1% to the emission observed by the 
LAT. We discuss these possibilities below. 

If the GeV and TeV 7-rays are not produced 
co-spatially (model assumption #1), this could ex- 
plain o ur results at high Bra. It is possible they 
are not iBottcher et al.|l2008l Essev et al.ll20lTl I2010L 
lEssev k Kusenkoll2010| . 



In this case one might ex- 
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Figure 3: Similar to Fig. d only for 1ES 1101-232. 
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Figure 4: Similar to Figs. [2]and[3l only here we show 
the combined results for both 1ES 1101—232 and 
1ES 0229+200. 



pect the TeV flux to be larger than the extrapolated 
LAT flux. If the two sources we used are variable (as- 
sumption #2), the variablity would have to be greater 
than the observational errors, which are quite large, 
to invalidate our results. Our results could be inter- 
preted as evidence for axion-like particles (assump- 
tion #3), which would allow TeV photons to avoid 
much of the EBL attenuation. This would be diffi- 
cult to distinguish, although with sufficiently accurate 
observations certain 7-ray signatures may be found 
[Sanchez-Conde et al 1 120091 e.g.,]. Finally, if all our 
other assumptions are correct, the cascade component 
is likely to exist at the > 5cx level. This implies that 
it is very unlikely that plasma beam instabilit ies would 
eliminate the cascade iBroderick et al.l 2012, as sump- 



tion #4] . See also IVenters fc Pavlidoul 2012j | for a 
critical assesment of this possibility. 

Of our model parameters, the EBL model of 
iFinke et all [201dj is unlikely to be wrong to a large 



degree, given its proximity to galaxy counts and agree- 
ment with other models. The timescale the blazar has 
been operating would not effect the upper limits on 
Big, nor could the jet opening angle. We will explore 
other EBL models in future work. 

The final possibility is that there is indeed a sig- 
nificant cascade component created, and the upper 
limits on Big an d Lb can be taken at face value. 
We note that in this case, the parameters are in the 
range where one might expect extended GeV emis- 
sion ("7-ray halos") around TeV sources to b e de- 
tectable by the LAT [Neronov fc Semikoj[2009| . The 
detection of this emission or lack thereof could pro- 
vide confirmation of these parameter constraints, or 
the other possibilities discussed above. The detection 
of t hese halos in the Fermi era has been controver- 
sial lAckermann et al. 1 12011 lAndo fc Kusenkdl2010L 
iNeronov et alJl201l| . 

The constraints on lower magnetic field values from 
the cascade are less certain, not just because of their 
lower significance but also due to their dependence 
on parameter assumptions. The timescale the blazar 
has been emitting TeV 7-rays or changes the opening 



angle can significant l y effe ct the cascade Arlen et al.l 
2012], iDermer et al.l |2?5TT |. We will explore how 



changes of these parameters effect our results. 

The IGMF may have been ge nerated in phase tran- 
sition s in the early universe [e.g. jGrasso fc Rubinstein] 
12001 1 . Our preliminary results ruling out low Big 
disfavor IGMFs generated from electro-weak phase 
transitions, if the genera ted fields are non-helical 
[Neronov fc Semikoz [2009]. However, as discussed 
above, the lower magnetic field values are dependent 
on caveats including those regarding source variabil- 
ity. If the IGMF is generated from phase transitions, 
and if it is possible to constrain Big > 10 -15 G, then 
it will be almost impo ssible to detect gravitationa l 
waves from inflation [Fujita fc Mukohvamal l2012j j. 
The IGMF parameters also have implications for the 
diffuse gamma ray background and its anisotrq py 
Ackermann et al. I l2012bl IVenters fc Pavlidou1l2012| . 
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